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Abstract. Dark matter has been introduced to explain many independent gravitational 

effects at different astronomical scales, in galaxies, groups of galaxies, clusters, superclustcrs 
and even across the full horizon. This review describes the accumulated astronomical, 
astrophysical, and cosmological evidence for dark matter. It is written at a non-specialist 
level and intended for an audience with little or only partial knowledge of astrophysics or 
cosmology. 
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1. Introduction 

According to the standard "Cold Dark Matter" (CDM) paradigm, dark matter originated 
from quantum fluctuations in the almost uniform, early Universe during a very short period 
of cosmic inflation, and emerged with negligible thermal velocities and a Gaussian and scale- 
free distribution of density fluctuations. It is one of the most important tasks of our time to 
establish whether this paradigm is supported by observations. 

In spite of the importance of this task, there is no review of the accumulated astronomical, 
astrophysical, and cosmological observational testimony on the CDM paradigm. The hterature 
summing such evidence usually refers only to the rotation curves of spiral galaxies, to the 
observations of strongly lensing galaxies and clusters, and to the determination of parameters 
in fits combining the Cosmic Microwave Background (CMB) with some other constraints. 
Quite recently the "Bullet cluster" gained fame due to the combination of lensing observations 
of its gravitational field and X-ray observations of the heated Intra-Cluster Medium (ICM). 
But there exists much more evidence. 

The purpose of this review is to summarize all the gravitational evidence of dark matter 
without entering into speculation on its true nature. Whether the correct explanation implies 
conventional matter, unconventional particles, or modifications to gravitational theory, all the 
observed gravitational effects described below have to be explained. Thus the use of the term 
"dark matter" is generic and does not imply that we are taking sides. 

Beginning from the first startling discoveries of Jan Hcndrik Oort in 1932 [1] of missing 
matter in the Galactic disk which has not been confirmed by modern observations (Sec. 2), 
and that of Fritz Zwicky in 1933 [2] of missing matter in the Coma cluster (Sec. 3), much 
later understood to be "dark matter" , we describe the kinematics of virially bound systems 
(Sec. 3) and rotating spiral galaxies (Sec. 4). Different testimonies of missing mass in groups 
and clusters are offered by the comparison of visible hght and X-rays (Sec. 5). The mass 
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to luminosity ratios of objects at all scales can be compared to what would be required to 
close the Universe (Sec. 6), and mass autocorrelation functions relate galaxy masses to dark 
halo masses (Sec. 7). An increasingly important method to determine the weights of galaxies, 
clusters and gravitational fields at large, independently of electromagnetic radiation, is lensing, 
strong as well as weak (Sec. 8). 

In radiation the most important tools are the temperature and polarization anisotropics 
in the Cosmic Microwave Background (Sec. 9), which give information on the mean density of 
both dark matter and baryonic matter as well as on the geometry of the Universe. The large 
scale structures of matter exhibit similar fluctuations with similar information in the Baryonic 
Acoustic Oscillations (BAO) (Sec. 10). The amplitude of the temperature variations in the 
CMB prove, that galaxies could not have formed in a purely baryonic Universe (Sec. 11). 
Simulations of large scale structures also show that dark matter must be present (Sec. 12). 
The best quantitative estimates of the density of dark matter come from overall parametric 
fits (to various cosmological models) of CMB data, BAO data, and redshifts of supernovae of 
type la (SNela) (Sec. 13). A particularly impressive testimony comes from combining lensing 
information on pairs of colliding clusters with X-ray maps (Sec. 14). We conclude (Sec. 15) 
with a few general comments on the nature of dark matter. 

2. Stars near the Galactic disk 

In 1932 the Dutch astronomer Jan Hendrik Oort analyzed the vertical motions of all known 
stars near the Galactic plane and used these data to calculate the acceleration of matter [1]. 
This amounts to treating the stars as members of a " star atmosphere" , a statistical ensemble 
in which the density of stars and their velocity dispersion defines a "temperature" from which 
one obtains the gravitational potential. This is analogous to how one obtains the gravitational 
potential of the Earth from a study of the atmosphere. The result contradicted grossly the 
expectations: the potential provided by the known stars was not sufficient to keep the stars 
bound to the Galactic disk, the Galaxy should rapidly be losing stars. Since the Galaxy 
appeared to be stable there had to be some missing matter near the Galactic plane, Oort 
thought, exerting gravitational attraction. This used to be counted as the first indication for 
the possible presence of dark matter in our Galaxy. 

The possibility that this missing matter would be non-baryonic could not even be thought 
of at that time. Note that the first neutral baryon, the neutron, was discovered by James 
Chadwick in the same year, in 1932. 

However, it is nowadays considered, that this does not prove the existence of dark matter 
in the disk. The potential in which the stars are moving is not only due to the disk, but rather 
to the totahty of matter in the Galaxy which is dominated by the Galactic halo. The advent of 
much more precise data in 1998 led Holmberg and Flynn [3] to conclude that no dark matter 
was present in the disk. 

Oort determined the mass of the Galaxy to be 10^^ Mgun, and thought that the 
nonluminous component was mainly gas. Still in 1969 he thought that intcrgalactic gas made 
up a large fraction of the mass of the universe [4]. The general recognition of the missing 
matter as a possibly new type of non-baryonic dark matter dates to the early eighties. 
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3. Virially bound systems 

Stars move in galaxies and galaxies in clusters along their orbits; the orbital velocities are 
balanced by the total gravity of the system, similar to the orbital velocities of planets moving 
around the Sun in its gravitational field. In the simplest dynamical framework one treats 
clusters of galaxies as statistically steady, spherical, self-gravitating systems of objects of 
average mass m and average orbital velocity v. The total kinetic energy E of such a system is 
then 

E = ]^Nmv'^ . (1) 
If the average separation is r, the potential energy of N{N — l)/2 pairings is 

C/ = -^iV(iV-i)^. (2) 

The virial theorem states that for such a system 

E = -U/2 . (3) 
The total dynamic mass M can then be estimated from and r from the cluster volume 

M = Nm^ -— . (4) 

An empirical formula often used in simulations of the density distribution of dark matter 
halos in clusters is 

where po is a normalization constant and < a < 3/2. Standard choices are a = 1 for the 
Navarro- Frenk- White profile [5] and the a = 3/2 profile of Moore & al. [6], both cusped at 
r = 0. 

However, Gentile & al. [7] have shown that cusped profiles are in clear conflict with data on 
spiral galaxies. Central densities are rather flat, scaling approximately as po 'T' luminous- "^^^ 
best-fit disk -|- NFW halo mass model fits the rotation curves poorly, it imphes an implausibly 
low stellar mass-to-light ratio and an unphysically high halo mass. 



3.1. The Coma cluster 

Historically, the second possible indication of dark matter, the first time in an object at a 
cosmological distance, was found by Fritz Zwicky in 1933 [2]. While measuring radial velocities 
of member galaxies in the Coma cluster (that contains some 1000 galaxies), and the cluster 
radius from the volume they occupy, Zwicky was the first to use the virial theorem to infer 
the existence of unseen matter. He was able to infer the average mass of galaxies within the 
cluster, and obtained a value about 160 times greater than expected from their luminosity (a 
value revised today), and proposed that most of the missing matter was dark. His suggestion 
was not taken seriously at first by the astronomical community, that Zwicky felt as hostile and 
prejudicial. Clearly, there was no candidate for the dark matter because gas radiating X-rays 
and dust radiating in the infrared could not yet be observed, and non-baryonic matter was 
unthinkable. Only some forty years later when studies of motions of stars within galaxies also 
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Figure 1. Density profile of matter components enclosed within a given radius r in the Coma 
cluster versus r/rmriai- From E. L. Lokas & G. A. Mamon [8]. 



implied the presence of a large halo of unseen matter extending beyond the visible stars, dark 
matter became a serious possibility. 

Zwicky found to his surprise that the orbital velocities were almost a factor of ten larger 
than expected from the summed mass of all galaxies belonging to the Coma cluster. He then 
concluded that in order to hold galaxies together the cluster must contain huge amounts of some 
non-luminous matter. Since that time modern observations have revised our understanding of 
the composition of clusters. Luminous stars represent a very small fraction of a cluster mass; 
in addition there is a baryonic, hot intracluster medium (ICM) visible in the X-ray spectrum. 
Rich clusters typically have more mass in hot gas than in stars; in the largest virial systems 
like the Coma the composition is about 85% dark matter, 14% ICM, and only 1% stars [8]. 

In modern applications of the virial theorem one also needs to model and parametrize 
the radial distributions of the ICM and the dark matter densities. In the outskirts of galaxy 
clusters the virial radius roughly separates bound galaxies from galaxies which may either be 
infalling or unbound. The virial radius r^r is conventionally defined as the radius within which 
the mean density is 200 times the background density. 

Matter accretion is in general quite well described within the approximation of the 
Spherical Collapse Model. According to this model, the velocity of the infall motion and 
the matter overdensity are related. Mass profile estimation is thus possible once the infall 
pattern of galaxies is known [9]. 

In Fig. 1 the Coma profile is fitted [8] with Eq. (5) with a = which describes a centrally 
finite profile which is almost fiat. The separation of different components in the core is not 
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well done with Eq. (5) because the Coma has a binary center like many other clusters [10]. 

The physical size of clusters is expected to depend on the mass, characterized by the mass 
concentration index Cg = r^ir/rs, roughly Cs oc M"'^'^. This can be tested by strong lensing of 
a sequence of halos. 
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Figure 2. Density profile of matter components in the cluster AC 114, enclosed within a 
given projected radius. From M. Sereno & at [14]. 

3.2. The AC 114 cluster 

Dark matter is usually dissected from baryons in lensing analyses by first fitting the lensing 
features to obtain a map of the total matter distribution and then subtracting the gas mass 
fraction as inferred from X-ray observations [11, 12]. The total mass map can then be 
obtained with parametric models in which the contribution from cluster-sized DM haloes can 
be considered together with the main galactic DM haloes [13]. Mass in stars and stellar 
remnants is estimated converting galaxy luminosity to mass assuming suitable stellar mass to 
light ratios. 

One may go one step further by exploiting a parametric model which has three kinds of 
components: cluster-sized dark matter haloes, galaxy-sized (dark plus stellar) matter haloes, 
and a cluster-sized gas distribution [10, 14]. The results of such an analysis of the dynamically 
active cluster AC 114 is shown in Fig. 2. 

3.3. The Local Group 

The Local Group is a very small virial system, dominated by two large galaxies, the M31 or 
Andromeda galaxy, and the Milky Way. The M31 exhibits blueshift, falling in towards us. 
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Evidently the Galaxy and M31, as well as all the minor galaxies in the Local Group, form a 
bound system which is oscillating. 

In this virial system the two large galaxies dominate the dynamics, so that it is not 
meaningful to define, as for the Coma, an average pairwise separation between galaxies, an 
average mass or an average orbital velocity. The total kinetic energy E is still given by the 
sum over all the group members, and the potential energy U by the sum over all the galaxy 
pairs, but here the pair formed by the M31 and the Milky Way dominates, and the pairings 
of the small members with each other are negligible. 

An interesting recent claim is, that the mass estimate of the Local Group is also affected 
by the accelerated expansion, the "dark energy". A. D. Chernin & al. [15] have shown that the 
the potential energy U is reduced in the force field of dark energy, so that the virial theorem 
for N masses rrii with baryocentric radius vectors rj takes the form 

E^-Uj + U2^ (6) 

where U is defined as in Eq. (3), and 

C^ = -^E".^? (7) 

is a correction which reduces the potential energy due to the background dark energy density 
p^. In the Local Group this correction to the mass appears to be quite substantial, of the order 
of 30% - 50%. 

The dynamical mass of the local group comes out to be 3.2 — 3.7 10^^ solar masses whereas 
the total visible mass of the Galaxy + M31 is only 2 10^^ solar masses. Thus there is a large 
amount of dark matter. 



4. Rotation curves of spiral galetxies 

Spiral galaxies are stable gravitationally bound systems in which visible matter is composed 
of stars and interstellar gas. Most of the observable matter is in a relatively thin disc, where 
stars and gas rotate around the galactic center on nearly circular orbits. If the circular velocity 
at radius r is v in a galaxy with mass M(r) inside r, the condition for stability is that the 
centrifugal acceleration v/r should equal the gravitational pull GM{r)/r'^, and the radial 
dependence of v would then be expected to follow Kepler's law 

-V^. («) 

The surprising result from mcasiucmcnts of galaxy-rotation curves is that the velocity does 
not follow this 1 / ^Jr law, but stays constant after attaining a maximum at about 5 kpc. The 
most obvious solution to this is that the galaxies are embedded in extensive, diffuse haloes of 
dark matter. If the enclosed mass inside the radius r, M(r), is proportional to r it follows that 
v{r) ~ constant. 

The rotation curve of most galaxies can be fitted by the superposition of contributions 
from the stellar and gaseous disks, the luminous components, both modeled by thin exponential 
disks, sometimes a bulge, and the dark halo, modeled by a quasi-isothermal sphere. The 
inner part is difficult to model because the density of stars is high, rendering observations of 
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Figure 3. Best disk-halo fits to the Universal Rotation Curve (dotted/dashed line: disc/halo). 
Each object is identified by the halo virial mass [16]. 

individual star velocities difficult. Thus the ffis are not unique, the relative contributions of 
disk and dark matter is model-dependent, and it is not even sure whether galactic disks do 
contain dark matter. Typically the dark matter is about half of the total mass. 

In Fig. 3 we show the rotation curves fitted for eleven well-measured galaxies [16]. One 
notes, that at very small radii the dark halo component is indeed much smaller than the 
luminous disk component. At large radii, however, the need for a halo of dark matter is 
obvious. On galactic scales, the contribution of dark matter generally dominates the total mass. 
Note the contribution of the baryonic component, negligible for small masses but increasingly 
important in the larger structures. 
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Our Galaxy is complicated because of what appears to be a density dip at 9 kpc and a 
small dip at 3 kpc, as is seen in Fig. 4 [17]. To fit the measured rotation curve one needs at 
least three contributing components: a central bulge, the star disk + gas, and a dark matter 
halo [17, 18, 19]. For small radii there is a choice of empirical rotation curves, and no dark 
matter component appears to be needed until radii beyond 15 kpc. 




R kpc 



Figure 4. Decomposition of the rotation curve of the Milky Way into the components bulge, 
stellar disk + interstellar gas, dark matter halo (the red curves from left to right). From Y. 
Sofue et al. [17]. 
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5. Small galaxy groups emitting X-rays 

There are examples of groups formed by a small mimbcr of galaxies which are enveloped in a 
large cloud of hot gas (ICM), visible by its X-ray emission. One may assume that the electron 
density distribution associated with the X-ray brightness is in hydrostatic equilibrium, and 
one can extract the ICM radial density profiles by fits. 

The amount of matter in the form of hot gas can be deduced from the intensity of this 
radiation. Adding the gas mass to the observed luminous matter, the total amount of baryonic 
matter, M^, can be estimated [20, 21]. In clusters studied, the gas fraction increases with the 
distance from the center; the dark matter is more concentrated than the visible matter. 

The temperature of the gas depends on the strength of the gravitational field, from which 
the total amount of gravitating matter, Mg^av, in the system can be deduced. In many such 
small galaxy groups one finds Mgrav/^b ^ 3. Thus a dark halo must be present. An accurate 
estimate of Mg^av requires that also dark energy is taken into account, because it reduces 
the strength of the gravitational potential. There are sometimes doubts whether all galaxies 
apparently in these groups are physical members. If not, they will artificially increase the 
velocity scatter and thus lead to larger virial masses. 

On the scale of large clusters of galaxies hke the Coma, it is generally observed that dark 
matter represents about 85% of the total mass and that the visible matter is mostly in the 
form of a hot ICM. 

6. Mass to luminosity ratios 

Let us define the mass/luminosity ratio of an astronomical object as T = M/L. Stellar 
populations exhibit values T = 1 — 10 in solar units, in the solar neighborhood T = 2.5 — 7, 
in the Galactic disk T = 1.0 — 1.7 [22]. In the Coma the value of T depends strongly on the 
halo profile for r < 0.2 Mpc; beyond that the value is 400 and slowly decreasing with radius 
[8]. Many dwarf spheroidal galaxies exhibit very high ratios: in Draco T = 330 ± 125, in 
Andromeda IX T = 93^^^°. The cluster AC 114 discussed in Sec. 2.2 is very underluminous 
having T = 700 ± 100 [14]. 

A special case is the ultra-faint dwarf disk galaxy Segue 1 [23] which has a baryon mass 
of only about 1000 solar masses. One interpretation is that it is a thin non-rotating stellar 
disk not accompanied by a gas disk, embedded in an axisymmetric DM halo and with an 
/ = Mhaio/Mi, fa 200. But if the disk rotates, / could be as high as 2000. If Segue 1 also has a 
magnetized gas disk, the DM halo has to confine the effective pressure in the stellar disk and 
the magnetic Lorentz force in the gas disk as well as possible rotation. Then / could be very 
large [23]. 

These mass/luminosity ratios are however local, characteristic for isolated objects but not 
for the whole observable Universe. The value for ^universe is roughly known. One can then 
deduce that the critical density to close the Universe would require T ^ 990. However, there 
is not sufficiently much luminous matter to achieve closure, therefore a large amount of DM is 
needed. 
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Figure 5. Dark matter halo mass Mhaio as a function of stellar mass . The thick black curve 
is the prediction from abundance matching assuming no dispersion in the relation between the 
two masses. Red and green dashed curves assume some dispersion in log A/*. The dashed 
black curve is the satellite fraction as a function of stellar mass, as labeled on the axis at the 
right-hand side of the plot. From Qi Guo & al. [24]. 

7. Mass autocorrelation functions 

If galaxy formation is a local process, then on large scales galaxies must trace mass. This 
requires the study of how galaxies populate dark matter haloes. In simulations one attempts 
to track galaxy and dark matter halo evolution across cosmic time in a physically consistent 
way, providing positions, velocities, star formation histories and other physical properties for 
the galaxy populations of interest. 

Guo & al. [24] use abundance matching arguments to derive an accurate relation between 
galaxy stellar mass and dark matter halo mass. They combine a stellar mass function 
based on spectroscopic observations with a precise halo/subhalo mass function obtained from 
simulations. Assuming this stellar mass halo mass relation to be unique and monotonic, they 
compare it with direct observational estimates of the mean mass of haloes surrounding galaxies 
of given stellar mass inferred from gravitational lensing and satellite galaxy dynamics data, 
and use it to populate haloes in simulations. The stellar mass halo mass relation is shown 
in Fig. 5. The implied spatial clustering of stellar mass turns out to be in remarkably good 
agreement with a direct and precise measurement. By comparing the galaxy autocorrelation 
function with the total mass autocorrelation function, as averaged over the Local Supercluster 
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Figure 6. Dark matter column density vs. dark matter halo mass in solar units. From A. 
Boyarsky & al. [25] 



(LSC) volume, one concludes that a large amount of matter in the LSC is dark. 

A similar study is that of Boyarsky & al. [25] who find a universal relation between DM 
column density and DM halo mass, satisfied by matter distributions at all observable scales, 
in halo sizes from 10* to 10^^ solar masses, as shown in Fig. 6. Such a universal property is 
difficult to explain without dark matter. 

8. Strong and weak lensing 

A consequence of the Strong Equivalence Principle (SEP) is that a photon in a gravitational 
field moves as if it possessed mass, and light rays therefore bend around gravitating masses. 
Thus celestial bodies can serve as gravitational lenses probing the gravitational field, whether 
baryonic or dark without distinction. 

Since photons are neither emitted nor absorbed in the process of gravitational light 
defiection, the surface brightness of lensed sources remains unchanged. Changing the size of 
the cross-section of a light bundle only changes the fiux observed from a source and magnifies it 
at fixed surface-brightness level. If the mass of the lensing object is very small, one will merely 
observe a magnification of the brightness of the lensed object an effect called microlensing. 
Microlensing of distant quasars by compact lensing objects (stars, planets) has also been 
observed and used for estimating the mass distribution of the lens-quasar systems. 

Weak Lensing refers to deflection through a small angle when the light ray can be treated 
as a straight line (Fig. 7), and the deflection as if it occurred discontinuously at the point of 



Dark Matter: The evidence from astronomy, astrophysics and cosmology 



13 




Figure 7. Wave fronts and light rays in the presence of a cluster perturbation. From N. 
Straumann [26]. 



closest approach (the thin- lens approximation in optics). One then only invokes SEP which 
accounts for the distortion of clock rates. 
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Figure 8. This image resulted from color-subtracteion of a lensing singular isothermal 
elliptical galaxy. The strongly lensed object forms two prominent arcs A, B and a less extended 
third image C. From R.J. Smith & al. [27] 



Dark Matter: The evidence from astronomy, astrophysics and cosmology 



14 




1H-B qUri imil IHg lUB Htm 



Figure 9. Map of the dark matter distribution in the 2-square degree COSMOS field [28]: the 
hnear blue scale on top shows the gravitational lensing magnification k, which is proportional 
to the projected mass along the line of sight. Contours begin at 0.4% and are spaced by 0.5% 
in K. 

In Strong Lensing the photons move along geodesies in a strong gravitational potential 
which distorts space as well as time, causing larger deflection angles and requiring the full 
theory of GR. The images in the observer plane can then become quite complicated because 
there may be more than one null geodesic connecting source and observer; it may not even be 
possible to flnd a unique mapping onto the source plane c.f. Fig. 7 . Strong lensing is a tool 
for testing the distribution of mass in the lens rather than purely a tool for testing GR. An 
illustration is seen in Fig. 8 where the lens is an elliptical galaxy [27]. 

At cosmological distances one may observe lensing by composed objects such as galaxy 
groups which are ensembles of "point-like" , individual galaxies. Lensing effects are very model- 
dependent, so to learn the true magniflcation effect one needs very detailed information on the 
structure of the lens. 

The large-scale distribution of matter in the Universe is inhomogeneous in every direction, 
so one can expect that everything we observe is displaced and distorted by weak lensing. Since 
the tidal gravitational fleld and the deflection angles depend neither on the nature of the matter 
nor on its physical state, light deflection probes the total projected mass distribution. Lensing 
in infrared light offers an additional advantage of being able to sense distant background 
galaxies, since their number density is higher than in the optical range. 

Background galaxies would be ideal tracers of distortions if they were intrinsically circular. 
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Figure 10. Strong lensing analysis of the galaxy cluster Abell 370 (z=0.375) shown with 
the locations of the identified multiple systems within the critical line at z—1.2 in white (the 
redshift for the majority of the mutiple images). The red line delimits the region of multiple 
images for very high redshift sources (here assuming z = 6). From J. Richard & al. [29] 



because lensing transforms circular sources into ellipses. Any measured ellipticity would then 
directly reflect the action of the gravitational tidal field of the interposed lensing matter, 
and the statistical properties of the distortions would reflect the properties of the matter 
distribution. But many galaxies are actually intrinsically elliptical, and the ellipses are 
randomly oriented. This introduces noise into the inference of the tidal field from observed 
ellipticities. A useful feature in the sky is a fine-grained pattern of faint and distant blue 
galaxies appearing as a 'wall paper'. This makes statistical weak-lensing studies possible, 
because it allows the detection of the coherent distortions imprinted by gravitational lensing 
on the images of the galaxy population. Thus weak lensing has become an important technique 
to map non-luminous matter. A reconstruction of the largest and most detailed weak lensing 
survey ever undertaken with the Hubble Space Telescope is shown in Fig. 9 [28]. This map 
covers a large enough area to see extended filamentary structures. 

In Fig. 10 we see a first detailed example of the strong lensing of a cluster [29]. 

9. Cosmic Microwave Background (CMB) 

The tight coupling between radiation and matter density before decoupling caused the 
primordial adiabatic perturbations to oscillate in phase. Beginning from the time of last 
scattering, the receding horizon has been revealing these frozen density perturbations, setting 
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Figure 11. The maps of CMB radiation temperature (TT) and temperatme-polarization 
correlations (TE) from WMAP show anisotropics which can be analyzed by power spectra as 
functions of multipole moments. From G. Hinshaw & at [30] 



up a pattern of standing acoustic waves in the baryon-photon fluid. After decoupling, this 
pattern is visible today as temperature anisotropics with a certain regularity across the sky. 

The primordial photons are polarized by the anisotropic Thomson scattering process, but 
as long as the photons continue to meet free electrons their polarization is washed out, and no 
net polarization is produced. At a photon's last scattering, however, the induced polarization 
remains and the subsequently free- streaming photon possesses a quadrupole moment. 

Temperature and polarization fluctuations are analyzed in terms of multipole components 
or powers. The resulting distribution of powers versus multipole i, or multipole moment 
k = 271 /£, is the power spectrum which exhibits conspicuous Doppler peaks. In Fig. 11 
we display the TT and TE power spectra from the five-year data of WMAP as functions 
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of multipole moments [30]. The spectra can then be compared to theory, and theoretical 
parameters determined. Many experiments have determined the power spectra, so a wealth of 
data exists. 




Figure 12. The E-mode polarization power spectrum (EE) from the CMB observations of 
the QUaD collaboration & al. [32] 



Baryonic matter feels attractive self-gravity and is pressure-supported, whereas dark 
matter only feels attractive self-gravity, but is pressureless. Thus the Doppler peaks in the 
CMBR power spectrum testify about baryonic and dark matter, whereas the troughs testify 
about rarefaction caused by the baryonic pressure. The position of the first peak determines 
^mh"^- Combining the 5-year WMAP measurements of the temperature power spectrum (TT) 
with determinations of the Hubble constant h, the WMAP team finds the total mass density 
parameter Qm ~ 0.26 [31]. The ratio of amplitudes of the second-to-first Doppler peaks 
determines the baryonic density parameter Qi, ^ 0.04; the dark matter component is then 
Qdm ~ 0.22. 

Information on the temperature - E-mode polarization correlations (TE) comes from 
several measurements, among them WMAP [30, 31], and on the E-mode polarization power 
spectrum alone (EE) from the QUAD collaboration [32], Fig. 12. 

The results show two surprises: Firstly, since Qm <^ 1, a large component Q\ ~ 0.74 is 
missing, of unknown nature, and termed dark energy. The second surprise is that ordinary 
baryonic matter is only a small fraction of the total matter budget. The remainder is then 
dark matter, of unknown composition. Of the 4.2% of baryons in the Universe only about 1% 
is stars. 



10. Baryonic acoustic oscillations (BAO) 

A cornerstone of cosmology is the Copernican principle, that matter in the Universe is 
distributed homogeneously, if only on the largest scales of superclusters separated by voids. On 
smaller scales we observe inhomogeneities in the forms of galaxies, galaxy groups, and clusters. 
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The common approach to this situation is to turn to non-relativistic hydrodynamics and treat 
matter in the Universe as an adiabatic, viscous, non-static fluid, in which random fluctuations 
around the mean density appear, manifested by compressions in some regions and rarefactions 
in other. The origin of these density fluctuations was the tight couphng estabhshed before 
decouphng between radiation and charged matter density, causing them to oscillate in phase. 
An ordinary fluid is dominated by the material pressure, but in the fluid of our Universe three 
effects are competing: gravitational attraction, density dilution due to the Hubble flow, and 
radiation pressure felt by charged particles only. 

The inflationary fluctuations crossed the post-inflationary Hubble radius, to come back 
into vision with a wavelength corresponding to the size of the Hubble radius at that moment. 
At time t^q the overdensities began to amplify and grow into larger inhomogeneities. In 
overdense regions where the gravitational forces dominate, matter contracts locally and 
attracts surrounding matter, becoming increasingly unstable until it eventually collapses into 
a gravitationally bound object. In regions where the pressure forces dominate, the fluctuations 
move with constant amplitude as sound waves in the fluid, transporting energy from one region 
of space to another. 

Inflationary models predict that the primordial mass density fluctuations should be 
adiabatic, Gaussian, and exhibit the same scale invariance as the CMB fluctuations. The 
baryonic acoustic oscillations can be treated similarly to CMB, they are specifled by the 
dimensionless mass autocorrelation function which is the Fourier transform of the power 
spectrum of a spherical harmonic expansion. The power spectrum is shown in Fig. 13 [33]. 

As the Universe approached decoupling, the photon mean free path increases and 
radiation can diffused from overdense regions into underdense ones, thereby smoothing out 
any inhomogeneities in the plasma. The situation changed dramatically at recombination, at 
time 380 000 yr after Big Bang, when all the free electrons suddenly disappeared, captured 
into atomic Bohr orbits, and the radiation pressure almost vanished. Now the baryon acoustic 
waves and the CMB continued to oscillate independently, but adiabatically, and the density 
perturbations which have entered the Hubble radius since then can grow with full vigor into 
baryonic structures. 

The scale of BAO depends on and on the Hubble constant, h, so one needs information 
on h to break the degeneracy. The result is then Qm ~ 0.26. In the ratio Qh/Qm the h- 
dependence cancels out, so one can also quantify the amount of dark matter on very large 
scales by flb/Vlrn = 0.18 ± 0.04. 

11. Galcixy formation in purely bciryonic matter 

We have seen in Sec. 9 that the baryonic density parameter, Qb, is very small. The critical 
density ^critical is determined by the expansion speed of the Universe, and the mean baryonic 
density of the Universe (stars, interstellar & intergalactic gas) is only VL^ = 0.046 x ^critical 
[31]. 

The question arises whether the galaxies could have formed from primordial density 
fluctuations in a purely baryonic medium. We have also noted, that the fluctuations in CMB 
and BAO maintain adiabaticity. The amplitude of the primordial baryon density fluctuations 
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Figure 13. BAO in power spectra calculated from (a) the combined SDSS and 2dFGRS 
main galaxies, (b) the SDSS DR5 LRG sample, and (c) the combination of these two samples 
(solid symbols with la errors) . The data are correlated and the errors arc calculated from the 
diagonal terms in the covariance matrix. A Standard ACDM distance-redshift relation was 
assumed to calculate the power spectra with Q^n = 0.25, = 0.75. From Percival & al. [S3]. 

would have needed to be very large in order to form the observed number of galaxies. But then 
the amphtude of the CMB fluctuations would also have been very large, leading to intolerably 
large CMB anisotropies today. Thus galaxy formation in purely baryonic matter is ruled out 
by this argument alone. 

Thus one concludes, that the galaxies could only have been formed in the presence of 
gravitating dark matter which started to fluctuate early, unhindered by radiation pressure. 
This conclusion is further strengthened in the next Section. 

12. Large Scale Structures simulated 

In the ACDM paradigm, the nonlinear growth of dark matter structure is a well-posed problem 
where both the initial conditions and the evolution equations are known (at least when the 
effects of the baryons can be neglected) . 

The Aquarius Project [34] is a Virgo Consortium programme to carry out high-resolution 
dark matter simulations of Milky- Way-sized halos in the ACDM cosmology. This project seeks 
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clues to the formation of galaxies and to the nature of the dark matter by designing strategies 
for exploring the formation of our Galaxy and its luminous and dark satellites. The galaxy 




Figure 14. The left panel shows the projected dark matter density at z = in a slice of 
thickness 13.7 Mpc through the full box (137 Mpc on a side) of the 900^ parent simulation. 
The right panel show this halo resimulated at a different numerical resolution. The image 
brightness is proportional to the logarithm of the squared dark matter density projected along 
the line-of-sight. The circles mark r^Q, the radius within which the mean density is 50 times 
the background density. From V. Springel & al. [34] 

population on scales from 50 kpc to the size of the observable Universe has been predicted by 
hierarchical ACDM scenarios, and compared directly with a wide array of observations. So far, 
the ACDM paradigm has passed these tests successfully, particularly those that consider the 
large-scale matter distribution and has led to the discovery of a universal internal structure for 
dark matter halos As was noted already in Sec. 11, the observed structure of galaxies, clusters 
and superclusters, as illustrated by Fig. 14, could not have formed in a baryonic medium 
devoid of dark matter. 

Given this success, it is important to test ACDM predictions also on smaller scales, not 
least because these are sensitive to the nature of the dark matter. Indeed, a number of serious 
challenges to the paradigm have emerged on the scale of individual galaxies and their central 
structure. In particular, the abundance of small dark matter subhalos predicted within CDM 
halos is much larger than the number of known satellite galaxies surrounding the Milky Way. 

13. Dark matter from overall fits 

In Sec. 9 we have seen that the WMAP-5 year CMB data together with the Hubble constant 
value testify about the existence of dark matter [30, 31]. In Sec. 10 we addressed the BAO 
data [33] with the same conclusion. In overall fits one combines these with supernova data 
(SN la) which offer a constraint nearly orthogonal to that of CMB in the — fim-plane. 
The Union compilation of 414 SN la, which reduces to 307 SNe after selection cuts, includes 
the recent large samples of SNe la from the Supernova Legacy Survey and ESSENCE Survey, 
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Figure 15. 68.3 %, 95.4 % and 99.7% confidence level contours on ftj^ and flm obtained from 
CMB, BAO and the Union SN set, as well as their combination (assuming w — -1). Note the 
straight line corresponding to a flat Universe with n/^ + n„-, = 1. From M. Kowalski & al. [35]. 
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the older data sets, as well as the recently extended data set of distant supernovae observed 
with HST. M. Kowalski & al. [35] present the latest results from this Union compilation and 
discuss the cosmological constraints and its combination with CMB and BAO measurements. 
The CMB constraint is close to the line Q\ + Qm = 1, whereas the supernova constraint is 
close to the line — 1-6 x Qrn = 0.2. The BAO data constrain Qrn, but hardly at all ^a- This 
is shown in Fig. 15. 

Defining the vacuum energy density parameter by 

Qk = l-^A-^m , (9) 

a flat Universe corresponds to Qk = 0. For a ACDM Universe with a cosmological constant 
responsible for dark energy, a simultaneous fit to the data sets gives 

= o.285i°:°?° ± 0.011, = -o.oog^^Kooa , (lo) 

where the first error is statistical and the second error systematic. This is the most precise 
determination at present of the total mass density parameter Qm- At the same time one notes 
that the Universe is consistent with being fiat. Subtracting Qh = 0.046 from = 0.285 one 
obtains the density parameter for dark matter, Qdm ~ 0.24. Assuming flatness, M. Kowalski 
& al. [35] find = 0.274 ± 0.016 ± 0.013. 
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Figure 16. The merging cluster 1E0657-558. On the right is the smaller 'bullet' cluster which 
has traversed the larger cluster. The colors indicate the X-ray temperature of the plasma: blue 
is coolest and white is hottest. The green contours are the weak lensing reconstruction of the 
gravitational potential of the cluster. From D. Clowe & al. [36] 

14. Merging galaxy clusters 

A direct empirical proof of the existence of dark matter is furnished by observations of 1E0657- 
558, a unique cluster merger [36]. Due to the collision of two clusters, the dissipationless stellar 
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component and the fluid-like X-ray emitting plasma are spatially segregated. The gravitational 
potential observed by weak and strong lensing does not trace the plasma distribution which is 
the dominant baryonic mass component, but rather approximately traces the distribution of 
galaxies, cf. Fig. 16. The center of the total mass is offset from the center of the baryonic mass 
peaks, proving that the majority of the matter in the system is unseen. In front of the smaller 
'bullet' cluster which has traversed the larger one, a bow shock is evident in the X-rays. 

Two other merging systems with similar characteristics have been seen, although with 
lower spatial resolution and less clear-cut cluster geometry. In Fig. 17 we show the post-merging 
galaxy cluster MACS J0025.4-1222 [11]. 




Figure 17. The F555W-F814W color composite of the cluster MACS J0025.4.1222. Overlaid 
in red contours is the surface mass density k from the combined weak and strong lensing mass 
reconstruction. The contour levels are linearly spaced with Ak = 0.1, starting at Ak = 0.5, 
for a fiducial source at a redshift of — >■ oo. The X-ray brightness contours (also linearly 
spaced) are overlaid in yellow and the I-band light is overlaid in white. The measured peak 
position and error bars for the total mass are given as a cyan cross [11]. 

15. Comments and conclusions 

What we have termed "dark matter" is generic for observed gravitational effects on all scales: 
galaxies, small and large galaxy groups, clusters and superclusters, CMB anisotropics over the 
full horizon, baryonic oscillations over large scales, and cosmic shear in the large-scale matter 
distribution. The correct explanation or nature of dark matter is not known, whether it implies 
conventional matter, unconventional particles, or modifications to gravitational theory, but 
gravitational effects prove its existence in some form. 
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Only a few per cent of the total mass of the Universe is found in stars and hydrogen 
clouds, and this amount of baryonic matter is well accounted for by nucleosynthesis. If there 
exist particles which were very slow at time t^q when galaxy formation started, they could be 
candidates for cold dark matter. If dark matter were composed of such new particles, they 
must have become non-relativistic much earlier than the leptons, and then decoupled from the 
hot plasma. 

Whenever laboratory searches discover a new particle, it must pass several tests in order 
to be considered a viable DM candidate: it must be neutral, compatible with constraints 
on self-interactions (essentially coUisionless) , consistent with Big Bang nucleosynthesis, and 
match the appropriate relic density. It must be consistent with direct DM searches and 
gamma-ray constraints, it must leave stellar evolution unchanged, and be compatible with 
other astrophysical bounds. 

Theories attempting to explain dark matter without new particles, with new interactions 
or modified gravity, likewise have the burden to explain all the observed gravitational effects 
described in here. Thus there remains much to be done. 
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1. Introduction 

According to the standard "Cold Dark Matter" (CDM) paradigm, dark matter originated 
from quantum fluctuations in the almost uniform, early Universe during a very short period 
of cosmic inflation, and emerged with negligible thermal velocities and a Gaussian and scale- 
free distribution of density fluctuations. It is one of the most important tasks of our time to 
establish whether this paradigm is supported by observations. 

In spite of the importance of this task, there is no review of the accumulated astronomical, 
astrophysical, and cosmological observational testimony on the CDM paradigm. The hterature 
summing such evidence usually refers only to the rotation curves of spiral galaxies, to the 
observations of strongly lensing galaxies and clusters, and to the determination of parameters 
in fits combining the Cosmic Microwave Background (CMB) with some other constraints. 
Quite recently the "Bullet cluster" gained fame due to the combination of lensing observations 
of its gravitational field and X-ray observations of the heated Intra-Cluster Medium (ICM). 
But there exists much more evidence. 

The purpose of this review is to summarize all the gravitational evidence of dark matter 
without entering into speculation on its true nature. Whether the correct explanation implies 
conventional matter, unconventional particles, or modifications to gravitational theory, all the 
observed gravitational effects described below have to be explained. Thus the use of the term 
"dark matter" is generic and does not imply that we are taking sides. 

Beginning from the first startling discoveries of Jan Hcndrik Oort in 1932 [1] of missing 
matter in the Galactic disk which has not been confirmed by modern observations (Sec. 2), 
and that of Fritz Zwicky in 1933 [2] of missing matter in the Coma cluster (Sec. 3), much 
later understood to be "dark matter" , we describe the kinematics of virially bound systems 
(Sec. 3) and rotating spiral galaxies (Sec. 4). Different testimonies of missing mass in groups 
and clusters are offered by the comparison of visible hght and X-rays (Sec. 5). The mass 
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to luminosity ratios of objects at all scales can be compared to what would be required to 
close the Universe (Sec. 6), and mass autocorrelation functions relate galaxy masses to dark 
halo masses (Sec. 7). An increasingly important method to determine the weights of galaxies, 
clusters and gravitational fields at large, independently of electromagnetic radiation, is lensing, 
strong as well as weak (Sec. 8). 

In radiation the most important tools are the temperature and polarization anisotropics 
in the Cosmic Microwave Background (Sec. 9), which give information on the mean density of 
both dark matter and baryonic matter as well as on the geometry of the Universe. The large 
scale structures of matter exhibit similar fluctuations with similar information in the Baryonic 
Acoustic Oscillations (BAO) (Sec. 10). The amplitude of the temperature variations in the 
CMB prove, that galaxies could not have formed in a purely baryonic Universe (Sec. 11). 
Simulations of large scale structures also show that dark matter must be present (Sec. 12). 
The best quantitative estimates of the density of dark matter come from overall parametric 
fits (to various cosmological models) of CMB data, BAO data, and redshifts of supernovae of 
type la (SNela) (Sec. 13). A particularly impressive testimony comes from combining lensing 
information on pairs of colliding clusters with X-ray maps (Sec. 14). We conclude (Sec. 15) 
with a few general comments on the nature of dark matter. 

2. Stars near the Galactic disk 

In 1932 the Dutch astronomer Jan Hendrik Oort analyzed the vertical motions of all known 
stars near the Galactic plane and used these data to calculate the acceleration of matter [1]. 
This amounts to treating the stars as members of a " star atmosphere" , a statistical ensemble 
in which the density of stars and their velocity dispersion defines a "temperature" from which 
one obtains the gravitational potential. This is analogous to how one obtains the gravitational 
potential of the Earth from a study of the atmosphere. The result contradicted grossly the 
expectations: the potential provided by the known stars was not sufficient to keep the stars 
bound to the Galactic disk, the Galaxy should rapidly be losing stars. Since the Galaxy 
appeared to be stable there had to be some missing matter near the Galactic plane, Oort 
thought, exerting gravitational attraction. This used to be counted as the first indication for 
the possible presence of dark matter in our Galaxy. 

The possibility that this missing matter would be non-baryonic could not even be thought 
of at that time. Note that the first neutral baryon, the neutron, was discovered by James 
Chadwick in the same year, in 1932. 

However, it is nowadays considered, that this does not prove the existence of dark matter 
in the disk. The potential in which the stars are moving is not only due to the disk, but rather 
to the totahty of matter in the Galaxy which is dominated by the Galactic halo. The advent of 
much more precise data in 1998 led Holmberg and Flynn [3] to conclude that no dark matter 
was present in the disk. 

Oort determined the mass of the Galaxy to be 10^^ Mgun, and thought that the 
nonluminous component was mainly gas. Still in 1969 he thought that intcrgalactic gas made 
up a large fraction of the mass of the universe [4]. The general recognition of the missing 
matter as a possibly new type of non-baryonic dark matter dates to the early eighties. 
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3. Virially bound systems 

Stars move in galaxies and galaxies in clusters along their orbits; the orbital velocities are 
balanced by the total gravity of the system, similar to the orbital velocities of planets moving 
around the Sun in its gravitational field. In the simplest dynamical framework one treats 
clusters of galaxies as statistically steady, spherical, self-gravitating systems of objects of 
average mass m and average orbital velocity v. The total kinetic energy E of such a system is 
then 

E = ]^Nmv'^ . (1) 
If the average separation is r, the potential energy of N{N — l)/2 pairings is 

C/ = -^iV(iV-i)^. (2) 

The virial theorem states that for such a system 

E = -U/2 . (3) 
The total dynamic mass M can then be estimated from and r from the cluster volume 

M = Nm^ -— . (4) 

A generally used empirical formula to describe the density distribution of dark matter 
halos in clusters is 

where po is a normalization constant and < a < 3/2. Standard choices are a = 1 for the 
Navarro- Frenk- White profile [5] and the a = 3/2 profile of Moore & al. [6], both cusped at 
r = 0. 



3.1. The Coma cluster 

Historically, the second possible indication of dark matter, the first time in an object at a 
cosmological distance, was found by Fritz Zwicky in 1933 [2]. While measuring radial velocities 
of member galaxies in the Coma cluster (that contains some 1000 galaxies), and the cluster 
radius from the volume they occupy, Zwicky was the first to use the virial theorem to infer 
the existence of unseen matter. He was able to infer the average mass of galaxies within the 
cluster, and obtained a value about 160 times greater than expected from their luminosity (a 
value revised today), and proposed that most of the missing matter was dark. His suggestion 
was not taken seriously at first by the astronomical community, that Zwicky felt as hostile and 
prejudicial. Clearly, there was no candidate for the dark matter because gas radiating X-rays 
and dust radiating in the infrared could not yet be observed, and non-baryonic matter was 
unthinkable. Only some forty years later when studies of motions of stars within galaxies also 
implied the presence of a large halo of unseen matter extending beyond the visible stars, dark 
matter became a serious possibility. 

Zwicky found to his surprise that the orbital velocities were almost a factor of ten larger 
than expected from the summed mass of all galaxies belonging to the Coma cluster. He then 
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Figure 1. Density profile of matter components enclosed within a given radius r in the Coma 
cluster versus r/rmriai- From E. L. Lokas & G. A. Mamon [7]. 



concluded that in order to hold galaxies together the cluster must contain huge amounts of some 
non-luminous matter. Since that time modern observations have revised our understanding of 
the composition of clusters. Luminous stars represent a very small fraction of a cluster mass; 
in addition there is a baryonic, hot intracluster medium (ICM) visible in the X-ray spectrum. 
Rich clusters typically have more mass in hot gas than in stars; in the largest virial systems 
like the Coma the composition is about 85% dark matter, 14% ICM, and only 1% stars [7]. 

In modern applications of the virial theorem one also needs to model and parametrize 
the radial distributions of the ICM and the dark matter densities. In the outskirts of galaxy 
clusters the virial radius roughly separates bound galaxies from galaxies which may either be 
infalling or unbound. The virial radius r^j^ is conventionally defined as the radius within which 
the mean density is 200 times the background density. 

Matter accretion is in general quite well described within the approximation of the 
Spherical Collapse Model. According to this model, the velocity of the infall motion and 
the matter overdensity are related. Mass profile estimation is thus possible once the infall 
pattern of galaxies is known [8]. 

In Fig. 1 the Coma profile is fitted [7] with Eq. (5) with a = which describes a centrally 
finite profile which is almost flat. The separation of different components in the core is not 
well done with Eq. (5) because the Coma has a binary center like many other clusters [9]. 

The physical size of clusters is expected to depend on the mass, characterized by the mass 
concentration index Cs = Tmr/'Tsi roughly Cs oc M~^'^. This can be tested by strong lensing of 
a sequence of halos. 
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Figure 2. Density profile of matter components in the cluster AC 114, enclosed within a 
given projected radius. From M. Sereno & al. [13]. 

3.2. The AC 114 cluster 

Dark matter is usually dissected from baryons in lensing analyses by first fitting the lensing 
features to obtain a map of the total matter distribution and then subtracting the gas mass 
fraction as inferred from X-ray observations [10, 11]. The total mass map can then be 
obtained with parametric models in which the contribution from cluster-sized DM haloes can 
be considered together with the main galactic DM haloes [12]. Mass in stars and stellar 
remnants is estimated converting galaxy luminosity to mass assuming suitable stellar mass to 
light ratios. 

One may go one step further by exploiting a parametric model which has three kinds of 
components: cluster-sized dark matter haloes, galaxy-sized (dark plus stellar) matter haloes, 
and a cluster-sized gas distribution [9, 13]. The results of such an analysis of the dynamically 
active cluster AC 114 is shown in Fig. 2. 

3.3. The Local Group 

The Local Group is a very small virial system, dominated by two large galaxies, the M31 or 
Andromeda galaxy, and the Milky Way. The M31 exhibits blueshift, falling in towards us. 
Evidently the Galaxy and M31, as well as all the minor galaxies in the Local Group, form a 
bound system which is oscillating. 

In this virial system the two large galaxies dominate the dynamics, so that it is not 
meaningful to define, as for the Coma, an average pairwise separation between galaxies, an 
average mass or an average orbital velocity. The total kinetic energy E is still given by the 
sum over all the group members, and the potential energy U by the sum over all the galaxy 
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pairs, but here the pair formed by the M31 and the Milky Way dominates, and the pairings 
of the small members with each other are negligible. 

An interesting recent claim is, that the mass estimate of the Local Group is also affected 
by the accelerated expansion, the "dark energy". A. D. Chernin & al. [14] have shown that the 
the potential energy U is reduced in the force field of dark energy, so that the virial theorem 
for masses mj with baryocentric radius vectors takes the form 

£; = -^c/ + c/2, (6) 

where U is defined as in Eq. (3), and 

U2^ ^ l^mri (7) 

is a correction which reduces the potential energy due to the background dark energy density 
p^. In the Local Group this correction to the mass appears to be quite substantial, of the order 
of 30% - 50%. 

The dynamical mass of the local group comes out to be 3.2 — 3.7 lO"*^^ solar masses whereas 
the total visible mass of the Galaxy + M31 is only 2 10^^ solar masses. Thus there is a large 
amount of dark matter. 



4. Rotation curves of spiral galfLxies 

Spiral galaxies are stable gravitationally bound systems in which visible matter is composed 
of stars and interstellar gas. Most of the observable matter is in a relatively thin disc, where 
stars and gas rotate around the galactic center on nearly circular orbits. If the circular velocity 
at radius r is u in a galaxy with mass M{r) inside r, the condition for stability is that the 
centrifugal acceleration v/r should equal the gravitational pull GM{r)/r'^, and the radial 
dependence of v would then be expected to follow Kepler's law 

(8, 

The surprising result from measurements of galaxy-rotation curves is that the velocity does 
not follow this l/^/r law, but stays constant after attaining a maximum at about 5 kpc. The 
most obvious solution to this is that the galaxies are embedded in extensive, diffuse haloes of 
dark matter. If the enclosed mass inside the radius r, M(r), is proportional to r it follows that 
v{r) w constant. 

The rotation curve of most galaxies can be fitted by the superposition of contributions 
from the stellar and gaseous disks, both modeled by thin exponential disks, sometimes a 
bulge, and the dark halo, modeled by a quasi-isothermal sphere. The inner part is difficult to 
model because the density of stars is high, rendering observations of individual star velocities 
difficult. Thus the fits are not unique, the relative contributions of disk and dark matter 
is model-dependent, and it is not even sure whether galactic disks do contain dark matter. 
Typically the dark matter is about half of the total mass. 

In Figs. 3 and 4 we show the rotation curves fitted for two well-measured galaxies, NGC 
6503 and NGC 3198 [15]. One notes, that at very small radn the dark halo component is 
indeed much smaller than the stellar disk component. At large radii, however, the need for a 
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halo of dark matter is obvious. On galactic scales, the contribution of dark matter generally 
dominates the total mass. 

Our Galaxy is complicated because of what appears to be a density dip at 9 kpc and a 
small dip at 3 kpc, as is seen in Fig. 5 [16]. To fit the measured rotation curve one needs at 
least three contributing components: a central bulge, the star disk + gas, and a dark matter 
halo [16, 17, 18]. For small radii there is a choice of empirical rotation curves, and no dark 
matter component appears to be needed until radii beyond 15 kpc. 
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Figure 3. Decomposition of the rotation curve of NGC3198 with a submaximal stellar disk. 
From D. Fuchs [15]. 
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Figure 4. Decomposition of the rotation curve of NGC6503 with a submaximal stellar disk. 
From D. Fuchs [15]. 
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Figure 5. Decomposition of the rotation curve of the Milky Way into the components bulge, 
stellar disk + interstellar gas, dark matter halo (the red curves from left to right). From Y. 
Sofue et al. [16]. 



Dark Matter: The evidence from astronomy, astrophysics and cosmology 



10 



5. Small galaxy groups emitting X-rays 

There are examples of groups formed by a small mimbcr of galaxies which are enveloped in a 
large cloud of hot gas (ICM), visible by its X-ray emission. One may assume that the electron 
density distribution associated with the X-ray brightness is in hydrostatic equilibrium, and 
one can extract the ICM radial density profiles by fits. 

The amount of matter in the form of hot gas can be deduced from the intensity of this 
radiation. Adding the gas mass to the observed luminous matter, the total amount of baryonic 
matter, M^, can be estimated [19, 20]. In clusters studied, the gas fraction increases with the 
distance from the center; the dark matter is more concentrated than the visible matter. 

The temperature of the gas depends on the strength of the gravitational field, from which 
the total amount of gravitating matter, Mg^av, in the system can be deduced. In many such 
small galaxy groups one finds Mgrav/^b ^ 3. Thus a dark halo must be present. An accurate 
estimate of Mg^av requires that also dark energy is taken into account, because it reduces 
the strength of the gravitational potential. There are sometimes doubts whether all galaxies 
apparently in these groups are physical members. If not, they will artificially increase the 
velocity scatter and thus lead to larger virial masses. 

On the scale of large clusters of galaxies hke the Coma, it is generally observed that dark 
matter represents about 85% of the total mass and that the visible matter is mostly in the 
form of a hot ICM. 

6. Mass to luminosity ratios 

Let us define the mass/luminosity ratio of an astronomical object as T = M/L. Stellar 
populations exhibit values T = 1 — 10 in solar units, in the solar neighborhood T = 2.5 — 7, 
in the Galactic disk T = 1.0 — 1.7 [21]. In the Coma the value of T depends strongly on the 
halo profile for r < 0.2 Mpc; beyond that the value is 400 and slowly decreasing with radius 
[7]. Many dwarf spheroidal galaxies exhibit very high ratios: in Draco T = 330 ± 125, in 
Andromeda IX T = 93^^^°. The cluster AC 114 discussed in Sec. 2.2 is very underluminous 
having T = 700 ± 100 [13]. 

A special case is the ultra-faint dwarf disk galaxy Segue 1 [22] which has a baryon mass 
of only about 1000 solar masses. One interpretation is that it is a thin non-rotating stellar 
disk not accompanied by a gas disk, embedded in an axisymmetric DM halo and with an 
/ = Mhaio/Mi, fa 200. But if the disk rotates, / could be as high as 2000. If Segue 1 also has a 
magnetized gas disk, the DM halo has to confine the effective pressure in the stellar disk and 
the magnetic Lorentz force in the gas disk as well as possible rotation. Then / could be very 
large [22]. 

These mass/luminosity ratios are however local, characteristic for isolated objects but not 
for the whole observable Universe. The value for ^universe is roughly known. One can then 
deduce that the critical density to close the Universe would require T ^ 990. However, there 
is not sufficiently much luminous matter to achieve closure, therefore a large amount of DM is 
needed. 
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Figure 6. Dark matter halo mass Mhaio as a function of stellar mass . The thick black curve 
is the prediction from abundance matching assuming no dispersion in the relation between the 
two masses. Red and green dashed curves assume some dispersion in log A/*. The dashed 
black curve is the satellite fraction as a function of stellar mass, as labeled on the axis at the 
right-hand side of the plot. From Qi Guo & al. [23]. 

7. Mass autocorrelation functions 

If galaxy formation is a local process, then on large scales galaxies must trace mass. This 
requires the study of how galaxies populate dark matter haloes. In simulations one attempts 
to track galaxy and dark matter halo evolution across cosmic time in a physically consistent 
way, providing positions, velocities, star formation histories and other physical properties for 
the galaxy populations of interest. 

Guo & al. [23] use abundance matching arguments to derive an accurate relation between 
galaxy stellar mass and dark matter halo mass. They combine a stellar mass function 
based on spectroscopic observations with a precise halo/subhalo mass function obtained from 
simulations. Assuming this stellar mass - halo mass relation to be unique and monotonic, they 
compare it with direct observational estimates of the mean mass of haloes surrounding galaxies 
of given stellar mass inferred from gravitational lensing and satellite galaxy dynamics data, 
and use it to populate haloes in simulations. The stellar mass - halo mass relation is shown 
in Fig. 6. The implied spatial clustering of stellar mass turns out to be in remarkably good 
agreement with a direct and precise measurement. By comparing the galaxy autocorrelation 
function with the total mass autocorrelation function, as averaged over the Local Supercluster 
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Figure 7. Dark matter column density vs. dark matter halo mass in solar units. From A. 
Boyarsky & al. [24] 



(LSC) volume, one concludes that a large amount of matter in the LSC is dark. 

A similar study is that of Boyarsky & al. [24] who find a universal relation between DM 
column density and DM halo mass, satisfied by matter distributions at all observable scales, 
in halo sizes from 10* to 10^^ solar masses, as shown in Fig. 7. Such a universal property is 
difficult to explain without dark matter. 

8. Strong and weak lensing 

A consequence of the Strong Equivalence Principle (SEP) is that a photon in a gravitational 
field moves as if it possessed mass, and light rays therefore bend around gravitating masses. 
Thus celestial bodies can serve as gravitational lenses probing the gravitational field, whether 
baryonic or dark without distinction. 

Since photons are neither emitted nor absorbed in the process of gravitational light 
defiection, the surface brightness of lensed sources remains unchanged. Changing the size of 
the cross-section of a light bundle only changes the fiux observed from a source and magnifies it 
at fixed surface-brightness level. If the mass of the lensing object is very small, one will merely 
observe a magnification of the brightness of the lensed object an effect called microlensing. 
Microlensing of distant quasars by compact lensing objects (stars, planets) has also been 
observed and used for estimating the mass distribution of the lens-quasar systems. 

Weak Lensing refers to deflection through a small angle when the light ray can be treated 
as a straight line (Fig. 8), and the deflection as if it occurred discontinuously at the point of 
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Figure 8. Wave fronts and light rays in the presence of a cluster perturbation. From N. 
Straumann [25]. 



closest approach (the thin- lens approximation in optics). One then only invokes SEP which 
accounts for the distortion of clock rates. 
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Figure 9. This image resulted from color-subtracteion of a lensing singular isothermal 
elliptical galaxy. The strongly lensed object forms two prominent arcs A, B and a less extended 
third image C. From R.J. Smith & al. [26] 
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Figure 10. Map of the dark matter distribution in the 2-square degree COSMOS field [27]: the 
hnear blue scale on top shows the gravitational lensing magnification k, which is proportional 
to the projected mass along the line of sight. Contours begin at 0.4% and are spaced by 0.5% 
in K. 

In Strong Lensing the photons move along geodesies in a strong gravitational potential 
which distorts space as well as time, causing larger deflection angles and requiring the full 
theory of GR. The images in the observer plane can then become quite complicated because 
there may be more than one null geodesic connecting source and observer; it may not even be 
possible to flnd a unique mapping onto the source plane c.f. Fig. 8 . Strong lensing is a tool 
for testing the distribution of mass in the lens rather than purely a tool for testing GR. An 
illustration is seen in Fig. 9 where the lens is an elliptical galaxy [26]. 

At cosmological distances one may observe lensing by composed objects such as galaxy 
groups which are ensembles of "point-like" , individual galaxies. Lensing effects are very model- 
dependent, so to learn the true magniflcation effect one needs very detailed information on the 
structure of the lens. 

The large-scale distribution of matter in the Universe is inhomogeneous in every direction, 
so one can expect that everything we observe is displaced and distorted by weak lensing. Since 
the tidal gravitational fleld and the deflection angles depend neither on the nature of the matter 
nor on its physical state, light deflection probes the total projected mass distribution. Lensing 
in infrared light offers an additional advantage of being able to sense distant background 
galaxies, since their number density is higher than in the optical range. 

Background galaxies would be ideal tracers of distortions if they were intrinsically circular. 
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Figure 11. Strong lensing analysis of the galaxy cluster Abell 370 (z=0.375) shown with 
the locations of the identified multiple systems within the critical line at z—1.2 in white (the 
redshift for the majority of the mutiple images). The red line delimits the region of multiple 
images for very high redshift sources (here assuming z = 6). From J. Richard & al. [28] 



because lensing transforms circular sources into ellipses. Any measured ellipticity would then 
directly reflect the action of the gravitational tidal field of the interposed lensing matter, 
and the statistical properties of the distortions would reflect the properties of the matter 
distribution. But many galaxies are actually intrinsically elliptical, and the ellipses are 
randomly oriented. This introduces noise into the inference of the tidal field from observed 
ellipticities. A useful feature in the sky is a fine-grained pattern of faint and distant blue 
galaxies appearing as a 'wall paper'. This makes statistical weak-lensing studies possible, 
because it allows the detection of the coherent distortions imprinted by gravitational lensing 
on the images of the galaxy population. Thus weak lensing has become an important technique 
to map non-luminous matter. A reconstruction of the largest and most detailed weak lensing 
survey ever undertaken with the Hubble Space Telescope is shown in Fig. 10 [27]. This map 
covers a large enough area to see extended filamentary structures. 

In Fig. 11 we see a first detailed example of the strong lensing of a cluster [28]. 

9. Cosmic Microwave Background (CMB) 

The tight coupling between radiation and matter density before decoupling caused the 
primordial adiabatic perturbations to oscillate in phase. Beginning from the time of last 
scattering, the receding horizon has been revealing these frozen density perturbations, setting 
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Figure 12. The maps of CMB radiation temperature (TT) and temperatme-polarization 
correlations (TE) from WMAP show anisotropics which can be analyzed by power spectra as 
functions of multipole moments. From G. Hinshaw & at [29] 



up a pattern of standing acoustic waves in the baryon-photon fluid. After decoupling, this 
pattern is visible today as temperature anisotropics with a certain regularity across the sky. 

The primordial photons are polarized by the anisotropic Thomson scattering process, but 
as long as the photons continue to meet free electrons their polarization is washed out, and no 
net polarization is produced. At a photon's last scattering, however, the induced polarization 
remains and the subsequently free- streaming photon possesses a quadrupole moment. 

Temperature and polarization fluctuations are analyzed in terms of multipole components 
or powers. The resulting distribution of powers versus multipole i, or multipole moment 
k = 27r/£, is the power spectrum which exhibits conspicuous Doppler peaks. In Fig. 12 
we display the TT and TE power spectra from the 5 year data of WMAP as functions 
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of multipole moments [29]. The spectra can then be compared to theory, and theoretical 
parameters determined. Many experiments have determined the power spectra, so a wealth of 
data exists. 




Figure 13. The E-mode polarization power spectrum (EE) from the CMB observations of 
the QUaD collaboration & al. [31] 



Baryonic matter feels attractive self-gravity and is pressure-supported, whereas dark 
matter only feels attractive self-gravity, but is pressureless. Thus the Doppler peaks in the 
CMBR power spectrum testify about baryonic and dark matter, whereas the troughs testify 
about rarefaction caused by the baryonic pressure. The position of the first peak determines 
^mh"^- Combining the 5-year WMAP measurements of the temperature power spectrum (TT) 
with determinations of the Hubble constant h, the WMAP team finds the total mass density 
parameter Qm ~ 0.26 [30]. The ratio of amplitudes of the second-to-first Doppler peaks 
determines the baryonic density parameter Qi, ^ 0.04; the dark matter component is then 
Qdm ~ 0.22. 

Information on the temperature - E-mode polarization correlations (TE) comes from 
several measurements, among them WMAP [29, 30], and on the E-mode polarization power 
spectrum alone (EE) from the QUAD collaboration [31], Fig. 13. 

The results show two surprises: Firstly, since Qm <^ 1, a large component Q\ ~ 0.74 is 
missing, of unknown nature, and termed dark energy. The second surprise is that ordinary 
baryonic matter is only a small fraction of the total matter budget. The remainder is then 
dark matter, of unknown composition. Of the 4.2% of baryons in the Universe only about 1% 
is stars. 



10. Baryonic acoustic oscillations (BAO) 

A cornerstone of cosmology is the Copernican principle, that matter in the Universe is 
distributed homogeneously, if only on the largest scales of superclusters separated by voids. On 
smaller scales we observe inhomogeneities in the forms of galaxies, galaxy groups, and clusters. 
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The common approach to this situation is to turn to non-relativistic hydrodynamics and treat 
matter in the Universe as an adiabatic, viscous, non-static fluid, in which random fluctuations 
around the mean density appear, manifested by compressions in some regions and rarefactions 
in other. The origin of these density fluctuations was the tight couphng estabhshed before 
decouphng between radiation and charged matter density, causing them to oscillate in phase. 
An ordinary fluid is dominated by the material pressure, but in the fluid of our Universe three 
effects are competing: gravitational attraction, density dilution due to the Hubble flow, and 
radiation pressure felt by charged particles only. 

The inflationary fluctuations crossed the post-inflationary Hubble radius, to come back 
into vision with a wavelength corresponding to the size of the Hubble radius at that moment. 
At time t^q the overdensities began to amplify and grow into larger inhomogeneities. In 
overdense regions where the gravitational forces dominate, matter contracts locally and 
attracts surrounding matter, becoming increasingly unstable until it eventually collapses into 
a gravitationally bound object. In regions where the pressure forces dominate, the fluctuations 
move with constant amplitude as sound waves in the fluid, transporting energy from one region 
of space to another. 

Inflationary models predict that the primordial mass density fluctuations should be 
adiabatic, Gaussian, and exhibit the same scale invariance as the CMB fluctuations. The 
baryonic acoustic oscillations can be treated similarly to CMB, they are specifled by the 
dimensionless mass autocorrelation function which is the Fourier transform of the power 
spectrum of a spherical harmonic expansion. The power spectrum is shown in Fig. 14 [32]. 

As the Universe approached decoupling, the photon mean free path increases and 
radiation can diffused from overdense regions into underdense ones, thereby smoothing out 
any inhomogeneities in the plasma. The situation changed dramatically at recombination, at 
time 380 000 yr after Big Bang, when all the free electrons suddenly disappeared, captured 
into atomic Bohr orbits, and the radiation pressure almost vanished. Now the baryon acoustic 
waves and the CMB continued to oscillate independently, but adiabatically, and the density 
perturbations which have entered the Hubble radius since then can grow with full vigor into 
baryonic structures. 

The scale of BAO depends on and on the Hubble constant, h, so one needs information 
on h to break the degeneracy. The result is then Qm ~ 0.26. In the ratio Qh/Qm the h- 
dependence cancels out, so one can also quantify the amount of dark matter on very large 
scales by flb/Vlrn = 0.18 ± 0.04. 

11. Galcixy formation in purely bciryonic matter 

We have seen in Sec. 9 that the baryonic density parameter, Qb, is very small. The critical 
density ^critical is determined by the expansion speed of the Universe, and the mean baryonic 
density of the Universe (stars, interstellar & intergalactic gas) is only VL^ = 0.046 x ^critical 
[30]. 

The question arises whether the galaxies could have formed from primordial density 
fluctuations in a purely baryonic medium. We have also noted, that the fluctuations in CMB 
and BAO maintain adiabaticity. The amplitude of the primordial baryon density fluctuations 
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Figure 14. BAO in power spectra calculated from (a) the combined SDSS and 2dFGRS 
main galaxies, (b) the SDSS DR5 LRG sample, and (c) the combination of these two samples 
(solid symbols with Itr errors). The data are correlated and the errors arc calculated from the 
diagonal terms in the covariance matrix. A Standard ACDM distance - redshift relation was 
assumed to calculate the power spectra with Q^n = 0.25, = 0.75. From Percival & al. [32]. 

would have needed to be very large in order to form the observed number of galaxies. But then 
the amphtude of the CMB fluctuations would also have been very large, leading to intolerably 
large CMB anisotropies today. Thus galaxy formation in purely baryonic matter is ruled out 
by this argument alone. 

Thus one concludes, that the galaxies could only have been formed in the presence of 
gravitating dark matter which started to fluctuate early, unhindered by radiation pressure. 
This conclusion is further strengthened in the next Section. 

12. Large Scale Structures simulated 

In the ACDM paradigm, the nonlinear growth of dark matter structure is a well-posed problem 
where both the initial conditions and the evolution equations are known (at least when the 
effects of the baryons can be neglected) . 

The Aquarius Project [33] is a Virgo Consortium programme to carry out high-resolution 
dark matter simulations of Milky- Way-sized halos in the ACDM cosmology. This project seeks 
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clues to the formation of galaxies and to the nature of the dark matter by designing strategies 
for exploring the formation of our Galaxy and its luminous and dark satellites. The galaxy 




Figure 15. The left panel shows the projected dark matter density at z = in a slice of 
thickness 13.7 Mpc through the full box (137 Mpc on a side) of the 900^ parent simulation. 
The right panel show this halo resimulated at a different numerical resolution. The image 
brightness is proportional to the logarithm of the squared dark matter density projected along 
the line-of-sight. The circles mark r^Q, the radius within which the mean density is 50 times 
the background density. From V. Springel & al. [33] 

population on scales from 50 kpc to the size of the observable Universe has been predicted by 
hierarchical ACDM scenarios, and compared directly with a wide array of observations. So far, 
the ACDM paradigm has passed these tests successfully, particularly those that consider the 
large-scale matter distribution and has led to the discovery of a universal internal structure for 
dark matter halos As was noted already in Sec. 11, the observed structure of galaxies, clusters 
and superclusters, as illustrated by Fig. 15, could not have formed in a baryonic medium 
devoid of dark matter. 

Given this success, it is important to test ACDM predictions also on smaller scales, not 
least because these are sensitive to the nature of the dark matter. Indeed, a number of serious 
challenges to the paradigm have emerged on the scale of individual galaxies and their central 
structure. In particular, the abundance of small dark matter subhalos predicted within CDM 
halos is much larger than the number of known satellite galaxies surrounding the Milky Way. 

13. Dark matter from overall fits 

In Sec. 9 we have seen that the WMAP-5 year CMB data together with the Hubble constant 
value testify about the existence of dark matter [29, 30]. In Sec. 10 we addressed the BAO 
data [32] with the same conclusion. In overall fits one combines these with supernova data 
(SN la) which offer a constraint nearly orthogonal to that of CMB in the — fim-plane. 
The Union compilation of 414 SN la, which reduces to 307 SNe after selection cuts, includes 
the recent large samples of SNe la from the Supernova Legacy Survey and ESSENCE Survey, 
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Figure 16. 68.3 %, 95.4 % and 99.7% confidence level contours on ftj^ and flm obtained from 
CMB, BAO and the Union SN set, as well as their combination (assuming w — -1). Note the 
straight line corresponding to a flat Universe with n/^ + n„-, = 1. From M. Kowalski & al. [34]. 
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the older data sets, as well as the recently extended data set of distant supernovae observed 
with HST. M. Kowalski & al. [34] present the latest results from this Union compilation and 
discuss the cosmological constraints and its combination with CMB and BAO measurements. 
The CMB constraint is close to the line Q\ + Qm = 1, whereas the supernova constraint is 
close to the line Qa — 1.6 x Qrn = 0.2. The BAO data constrain Qrn, but hardly at all ^a- This 
is shown in Fig. 16. 

Defining the vacuum energy density parameter by 

Qk = l-^A-^m , (9) 

a flat Universe corresponds to Qk = 0. For a ACDM Universe with a cosmological constant 
responsible for dark energy, a simultaneous fit to the data sets gives 

= o.285i°:°?° ± 0.011, = -o.oog^^Kooa , (lo) 

where the first error is statistical and the second error systematic. This is the most precise 
determination at present of the total mass density parameter Qm- At the same time one notes 
that the Universe is consistent with being fiat. Subtracting Qh = 0.046 from = 0.285 one 
obtains the density parameter for dark matter, Qdm ~ 0.24. Assuming flatness, M. Kowalski 
& al. [34] find = 0.274 ± 0.016 ± 0.013. 
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Figure 17. The merging cluster 1E0657-558. On the right is the smaller 'bullet' cluster which 
has traversed the larger cluster. The colors indicate the X-ray temperature of the plasma: blue 
is coolest and white is hottest. The green contours are the weak lensing reconstruction of the 
gravitational potential of the cluster. From D. Clowe & al. [35] 

14. Merging galaxy clusters 

A direct empirical proof of the existence of dark matter is furnished by observations of 1E0657- 
558, a unique cluster merger [35]. Due to the collision of two clusters, the dissipationless stellar 
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component and the fluid-like X-ray emitting plasma are spatially segregated. The gravitational 
potential observed by weak and strong lensing does not trace the plasma distribution which is 
the dominant baryonic mass component, but rather approximately traces the distribution of 
galaxies, cf. Fig. 17. The center of the total mass is offset from the center of the baryonic mass 
peaks, proving that the majority of the matter in the system is unseen. In front of the smaller 
'bullet' cluster which has traversed the larger one, a bow shock is evident in the X-rays. 

Two other merging systems with similar characteristics have been seen, although with 
lower spatial resolution and less clear-cut cluster geometry. In Fig. 18 we show the post-merging 
galaxy cluster MACS J0025.4-1222 [10]. 




Figure 18. The F555W-F814W color composite of the cluster MACS J0025.4.1222. Overlaid 
in red contours is the surface mass density k from the combined weak and strong lensing mass 
reconstruction. The contour levels are linearly spaced with Ak = 0.1, starting at Ak = 0.5, 
for a fiducial source at a redshift of — >■ oo. The X-ray brightness contours (also linearly 
spaced) are overlaid in yellow and the I-band light is overlaid in white. The measured peak 
position and error bars for the total mass are given as a cyan cross [10]. 

15. Comments and conclusions 

What we have termed "dark matter" is generic for observed gravitational effects on all scales: 
galaxies, small and large galaxy groups, clusters and superclusters, CMB anisotropics over the 
full horizon, baryonic oscillations over large scales, and cosmic shear in the large-scale matter 
distribution. The correct explanation or nature of dark matter is not known, whether it implies 
conventional matter, unconventional particles, or modifications to gravitational theory, but 
gravitational effects prove its existence in some form. 
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Only a few per cent of the total mass of the Universe is found in stars and hydrogen 
clouds, and this amount of baryonic matter is well accounted for by nucleosynthesis. If there 
exist particles which were very slow at time t^q when galaxy formation started, they could be 
candidates for cold dark matter. If dark matter were composed of such new particles, they 
must have become non-relativistic much earlier than the leptons, and then decoupled from the 
hot plasma. 

Whenever laboratory searches discover a new particle, it must pass several tests in order 
to be considered a viable DM candidate: it must be neutral, compatible with constraints 
on self-interactions (essentially coUisionless) , consistent with Big Bang nucleosynthesis, and 
match the appropriate relic density. It must be consistent with direct DM searches and 
gamma-ray constraints, it must leave stellar evolution unchanged, and be compatible with 
other astrophysical bounds. 

Theories attempting to explain dark matter without new particles, with new interactions 
or modified gravity, likewise have the burden to explain all the observed gravitational effects 
described in here. Thus there remains much to be done. 
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